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Un  Diffraction  by  a Strip 
A.  Erdelyi  and  C.  H.  Papas 

Introduction 

Hie  problem  of  diffraction  by  an  infinite  strip  or  slit  has  been  the 

subject  of  several  investigations^.  There  are  at  least  two  "exact"  methods 

(2) 

for  attacking  this  problem.  One  of  these  is  the  integral  equation  method,'  ' 
the  other  the  Fourier-Lame  method. ^ The  integral  equation  obtained  for 
chis  problem  cannot  be  solved  in  closed  formj  expansion  of  the  solution  in 
powers  of  the  ratio  (strip  width/wavelength)  leads  to  useful  formulas  for  low 
frequencies.  In  the  Fourier-Lame  method  the  wave  equation  is  separated  in 
coordinates  of  the  elliptic  cylinder,  the  solution  appears  as  an  infinite 
scries  of  Mathieu  functions,  and  the  usefulness  of  the  result  is  limited  by 
the  convergence  of  these  infinite  series,  and  by  the  available  tabulation  of 
mathieu  functions. 

The  variational  technique  developed  by  Levine  and  Schwinger  avoids 
some  of  the  difficulties  of  the  above-mentioned  methods  and,  at  least  in 
principle,  is  capable  of  furnishing  good  approximations  for  all  frequency- 
ranges.  The  scattered  field  may  be  represented  as  the  effect  of  the  current 
induced  in  the  strip,  and  it  has  been  proved  by  Levine  and  Schwinger  ^ that 
it  is  possible  to  represent  the  amplitude  of  the  far-zone  scattered  field  in 
terms  of  the  induced  current  in  a form  which  is  stationary  with  respect  to 
small  variations  of  the  current  about  the  true  current.  Substitution,  in 
this  representation,  of  a rough  approximation  for  the  current  may  give  a re- 
markably good  approximation  of  the  far-zone  scattered  field  amplitude.  In 
this  note  we  assume  a normally  incident  field  polarized  parallel  to  the  gen- 
erators of  the  strip.  As  a rough  approximation,  we  take  a Uniterm  density  of 

the  current  induced  in  the  strip.  Since  the  incident  magnetic  field  is 

(?) 

constant  over  the  strip.  Fock’s  theory'^'  may  be  cited  in  support  of  the 


uniforraity  of  the  current  distribution,  except  near  the  edges  where  the  be- 
haviour of  the  field '■  indicates  an  infinite  current  density*  A mors 


(3) 

detailed  analysis  ol'  tne  current,  by  Houlrin  and  rhj.ilj.pa  , Is  available 
but  was  not  used  here . 

Once  the  (approximate)  amplitude  of  the  far-zone  field  has  been  ob- 
tained, the  scattering  cross-section  may  be  found  by  the  application  of  the 
scattering  theorem (k»7»R)  whiCfo  rQiates  this  cross-section  to  the  imaginary 
part  of  the  amplitude  of  the  far-zone  scattered  field  along  the  central  line 
of  the  umbral  region.  In  spite  of  the  crude  approximation  adopted  for  the 
induced  current,  the  scattering  cross-section  shows  a fair  agreement  with 
other  available  results. 


Integral  Equation 

We  assume  a plane  wave  with  harmonic  time  dependence  exp(-iut), 
normally  incident  on  the  perfectly  conducting  infinitesimally  thin  strip 
of  width  2a, 

z ■ 0,  -a^y^a  . (l) 

We  further  assume  that  the  incident  wave  is  polarized  parallel  to  the  edges 
of  the  strip  (i.e.  to  the  x axis)  so  that  the  only  nonvanishing  components 
of  its  electric  field  E,  magnetic  field  H,  and  complex  Poynting  vector 
(all  measured  in  KKS  units,  bars  denoting  conjugate  complex  numbers)  are 

Ex*  - exp(ikz)  (2) 

KyT  ■ k/(»p)  exp(ikz)  (3) 

sz  " Ejc*  ^ “ k/(«l*)  (U) 

The  scattered  electric  field  is  again  parallel  to  the  x axis,  and  the 
total  electric  field  is  the  sum  of  the  incident  and  scattered  field, 


E^Cy.z)  ♦ Exsc(y,z) 


(5) 


The  boundary  condition  on 


screen  is  the  vanishing  of  the  electric  field 


the 


^(y,©)  = 0 for 


'JC 


t<y 


(6) 


and  the  scattered  field,  Ej.  , must  represent,  at  large  distances  from  the 
screen,  an  outgoing  cylindrical  wave  (Sonmerfeld's  radiation  condition). 

The  expression  of  the  total  electric  field  in  terms  of  the  induced 
current  K^y)  is 

n 

|l/2. 


r (1)  r o 9"j1/2 

Ex(y*z)  ■ exp(ikz)  - o*i/U  J Hq'  (k  [(y  - y')  ♦ zdj  ) K^y'Jdy' 


-a 

ll/2. 


r?) 


) r 2 ^ 

i/k  H v (k  [(y-y')  + z J ) being  the  two-dimensional  free-space  Green's 

function.  The  scattered  field  in  (?)  certainly  satisfies  Sommerfeld 's  radia- 
tion condition!  in  order  that  it  also  satisfies  (6),  we  must  have 

(8) 


1 - cop/u  j"  (k  | y -y*  | ) K^y')  dy*  |y|^a 


and  this  is  the  integral  equation  of  our  problem. 


Far-Zone  Scattered  Field  Amplitude 
W9  define  the  far-zone  scattered  field  amplitude  A(0)  by 

Exsc^-i/U  [2/(n  i k p )]  exp(ikp)  A(0)  kp > (9) 

where  pcos  0»z  and  p sin  0 *■ -yj  thi3  expression  represents  a cylindrical 
outgoing  wave  of  "amplitude"  A. 

Since 

H ^ (k[(y-y' )2  ♦ z2]  ) fe/fnikp)!  ^ exp(ikp  ♦ iky ' sin  0) 

O \ *•  J / W J ftn^ 

when  kp — >00,  we  have  from  (7),  for  large  kp, 

Exsc— i/li  [2/(nikp|  exp(ikp)  J i«n  exp(iky'  sin  0)  \(y')  dy‘  (H) 


Comparing  (9)  and  (11)  we  get 

5 

A(p)  « ia>p  | exp(iky'  sin  0)  K_(v’)  (12) 

-a 

A Scattering  Theorem 

If  ? denotes  the  scattered  energy  flux  per  unit  length,  then  the 
scattering  cross-section  <r'/  is  defined  by 


cr'  - P/s  (13) 

where  S « 1/2  Re  S*-  k/(2»p)  is  the  incident  energy  flux  per  unit  area. 


a 


P - 1/2  Re  j"  (y!,°)  Kx(y’)  dy*  » 1/2  Re  J K^y')  d ly*  (1U) 

-a  -a 

since  E,/  =1  at  a ■ 0.  Hence, 

(15) 

-a 

But  from  (12) 


cr*  * top  A R0  ^ K (y’)  dy‘ 
U X 

-a 


A(0)  - i<un  J Kx(y»)  dy ! 

-a 

Therefore,  comparing  (15 ) and  (16)  we  get 


(16) 


CT1  * Im  A(0)/k 


(if) 


It  may  be  remarked  that  in  the  limit  for  very  high  frequencies,  the 
cross-section  (17)  turns  out  to  be  twice  the  projected  area.  At  first  sight 
the  factor  two  may  seem  inconsistent  with  geometrical  optics,  but  actually 
the  discrepancy  is  due  to  the  different  definitions  of  the  scattering  cross 
section. ^ 

Variational  Principle 

We  multiply  the  integral  equation  (8)  by  K^(y),  integrate  from  y*-a 


-5- 


tc  y-'ia,  divide  by  the  square  of  K^C>  ,‘dy,  and  then  recall  (16),  Thus 
v»s  cbtrj  n 


A(0)  - hi 


’ a 

J* Kx(y)  dy 

L-a 


a a 


J~  J Kx(y)  Ho(l)(k  jy-y'J)  yy«)  dy  dy' 


(18) 


-a  -a 


In  Appendix  A it  is  shown  that  the  expression  (18)  is  stationary  with  res- 
pect to  small  variations  of  about  the  true  K which  satisfies  the  integral 
equation  (8).  Consequently,  if  we  substitute  a reasonable  trial  function  for 
the  unknown  K,^  in  (18),  we  expect  to  obtain  a reasonably  accurate  value  of 
A(0). 


Trial  Function 

According  to  Fock's  theory'"'  the  current  induced  in  the  central 
region  of  the  strip  is  equal  to  twice  the  tangential  component  of  the  in- 
cident magnetic  field.  And  according  to  Bouwkamp^  Hy  must  have  a singu- 
larity like  (a^-y2 )”*/*■  at  the  edges.  This  certainly  conforms  to  the 

(3) 

detailed  information  Moullin  and  Phillips  reported.  In  the  present  paper 
we  take  for  a trial  function  K^y)  ■ 1,  thus  ignoring  the  effect  of  the 
edges . 


Scattering  Cross  Section 

With  K (v)  “ 1,  computation  (see  Appendix  B)  shows  that  (18)  be- 


comes 


f. 


2ka 


kG)!*1  . (Ui  a k)”  | f (t)dt-H.(x)  (2ka)  + (trike  T3/  (19) 


J 
V.  o 


And  applying  the  theorem  (17)  to  (19),  we  get 


-£>- 


?ka 


cr^/(Ua)  “ 


J .Jo  (t)di  - J3(?ka) 
o 


(w/ 


-2ka 

£ 

'2ka 

J JQ  (t)dt-J1(2ka) 
. 0 

♦ 

f I0(t)dt-yi(2ka)-(nkar1 

,5o 

The  integrals 


apy^aring  here  have  been  tabulated 


.....  do) 


Results 

A plot  of.  (T'/(hu)  versus  ka  is  shown  in  the  graph.  This  is  a 
plot  of  eq.(20)«  For  small  ka  we  have  Rayleigh  scattering.  As  ka  in- 
creases the  curve  performs  a damped  oscillation  about  <rV(U a)  * 1.  And 
for  ka  — ►°<’it  can  be  shown  by  means  of  asymptotic  representations  of  Bessel 
functions  and  their  integrals  that  a*— *l*a. 

IP  \ 

According  to  Babinet's  principle v/  the  problem  we  have  dis- 
cussed is  complementary  to  the  problem  of  scattering  by  a slit  for  a nor- 
mally incident  wavs  polarized  perpendicular  to  the  axis  of  the  slit.  For 
intermediate  values  of  ka  our  curve  behaves  qualitatively  as  Morse  and 
Rub  ins  tain's  the  quantitative  agreement  is  not  very  good,  the  devia- 

tion being  due  to  cur  choice  of  an  overly  simplified  trial  function.  It 
is  quite  remarkable  how  such  a rough  approximation  of  the  induced  current 
yields  fairly  good  results  over  the  entire  spectrum. 
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Appendix  A 


We  show  uere  that  A(0)  of  (18)  is  stationary  with  respect  tc  Small 
variations  of  Kx(y)  about  the  true  which  satisfies  the  integral 

equation  (8). 

From  (18)  and  the  calculus  of  variations  it  follows  that 
a a 

6 A f I"  yy)  «0(1)  (k  ly-y'i)  y.y)  ay  dy» 


J 

-a  -a 


* A ( f 5 My)  *V  (k  ly-y’l)  K,(y«)  dy  dy' 
•J  J 
-a  -a 


a a 


♦ A f J yy)  H0^j)  (k  Jy-y‘1)  6 Kx(y')  dy  dy»  - 


-a  -a 


a a 

2 J 6 yy)  dy  J yy)  dy  . 


-a 


-a 


fie  arranging  terms  we  get 


a 

6 A J Kx(y)  H0^  (k  I y-y'l)  Kx(y')  dy  dy’  - 

-a 

J 6 yy)  dy  j f yy)  dy  - A J H^1)  (k  | y-y’|)  yy')  dy'  ^ 

-a  -a  -a  J 

a f V V / 

)dy  * j j Kx(y)dy  - A J H^1)  (k|y-y'|)  Kx(y)  dy  (-  . 

(.-a  -a  ' 


* f 6 \()r’ 


The  right  side  of  this  equation  disappears  for  small  arbitrary  variations 

6 yy)  because  the  quantities  In  curly  brackets  are  zero  by  virtue  of  (12), 
x a 

i.e*,  A(0)  - lo>|i  \ yy)dy,  and  the  integral  equation  (8).  Consequently 

xj 

-A 

6 A ■ 0 and  the  expression  (18)  is  stationary. 


Appendix  E 


For  trial  function  K^(y)  - constant,  (io)  becomes 


wr 


/■  /* 

I6il2  J J H'^(1)  (k!y"y‘i)dy  dy' 

-a  -a 


Put  k(y  ♦ y * ) “S,  k (y-y»)  * t . Then 

*foT'  5TTK)2  [jHo(1)  (!t|> 


ds  dt 


where  S is  th«  square  with  vertices  (♦  2ka,0),  (0,^  2ka) 
Tot  Q be  the  quadrant  s >0,  t>0,  s ♦ t^2ka.  Then 


A(0)  8i(ka)2 


H U)  (t)  dt 


ilKB. 

(ka)2  f 


(2ka-t)H0(1)(t)dt. 


C < 2i 

Since  ^ t K0^(t)dt  - 2ka  (2ka) +~  , the  above  equation  becomes 


m3)  " UlkT^f  - Hi(1)(2K*>  * Tit-] 


J • 


And  this  is  (19). 


♦ 2ka 


*»  x 

a\ 


- 2ka 


♦ 2ka 


y 8 


-2ka 
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